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ABSTRACT 

We study the reliability of the statistical background subtraction method for com- 
puting the Ks-band luminosity function of cluster galaxies at z ~ 1 using mock 
Red-sequence Cluster Survey cluster catalogues constructed from GALFORM semi- 
analytic galaxies. The underlying cluster luminosity function in the mocks are compat- 
ible with recent estimates at z ^ 1 by several authors. We simulate different samples 
where the number of clusters with Ks-band photometry goes from 5 to a maximum 
of 50, in order to find the most suitable observational sample to carry out this study; 
the current observational status in the nIR wavelength range has been reached using 
5 real clusters at z ~ 1. We compute the composite luminosity function for several 
samples of galaxy clusters with masses ~ 1.5 x lO^^M© assuming a flux limited, com- 
plete sample of galaxies down to Ks = 21.0 magnitudes. We find that the Schechter fit 
parameters Ks* and a for a sample of galaxies with no redshift information are rather 
poorly constrained if both parameters are allowed to vary freely; if a is fixed at a 
fiducial value, then Ks* shows significantly improved stochastic uncertainties but can 
be influenced by systematic deviations. We find a significantly improved accuracy in 
the luminosity function parameters when adding photometric redshift information for 
bright cluster galaxies. The impact of a ten-fold increase in the number of clusters with 
available Ks-band photometry is that of decreasing stochastic errors in Ks* and a by 
factors of ~ 2 and ~ 4, respectively, for accuracies of up to (AKs*/Ks*)g^Q^j^^gjj^ ~ 0.02 
and (AQ!/Q;)gjQj.jjjjgjip — 0.09. The dwarf-to-giant ratios inferred from the luminosity 
functions of red-sequence galaxies in the mock catalogue agree very well with the 
underlying values; however, there is an indication that the semi-analytic model over- 
predicts the abundance of dwarf galaxies by up to a factor of 3 with respect to recent 
measurements. Finally, we find that in order to use estimates of Ks* to study the 
formation redshift of cluster galaxies at z = 1, the sample would need to contain 520 
z ~ 1 clusters, for an accuracy of ^ 2 Gyr at the 68 per cent confidence level. However, 
combining this method with other estimates may reduce significantly the sample size, 
and allow important new constraints on galaxy formation models. 
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1 INTRODUCTION 

The galaxy luminosity function (LP) is one of the fundamen- 
tal quantities of observational cosmology. Its evolution with 
redshift and dependence on galaxy morphology and environ- 
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ment has been extensively used to provide strong constraints 
on galaxy evolution models. 

Clusters of galaxies are the largest virialised structures 
at any given epoch, which grow from high-density regions 
in the primordial matter distribution of the Universe, and 
may have formed be tween 9 Gyrs ago and the present 
I Stanford et al.l l2006l 'l . Most of them have been detected 
in surveys carried out in the X-ray IjRomer et al. I I2OOII , 
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Barkhouse et al. 20061) and optical bands (Postman et alj 
19961 . iGonzalez et al.M200ll . iGladders fc Yeel |2005| ). and it 



is expected that upcoming Sunyaev-Zel' dovich surveys wi ll 
yield thousands of galaxy clusters (Carlstrom et al.]|2002l ). 
The fact that clusters contain large numbers of galaxies 
practically at the same distance from the observer within 
a small area on the sky, makes them ideal systems to study 
the LF down to very faint magnitudes. 

In the local Universe, the LF of clus t er ga laxies 
has been studie d exte nsively by iGoto et al.l l|2002t) and 
iDe Propris et al.1 (|2003l ). who found that there is no sig- 
nificant evidence for variations in the LF across a broad 
range of cluster properties, and also, that the cluster LF 
differs significantly from the field LF. Although many ef- 
forts have been done to perform such a study at z ~ 1, 
a significant look-back time of 7 Gyrs in the currently 
favoured cosmology, only a handful of high redshift clus- 
ters have been detected to date in optical and X-ray sur- 



veys (jStanford et al.lll997l.lToft et al.ll20q4lDe Propris et all 



I2OO7I ) . Recently. lEisenhardt et al.l l|2008l ) found 106 galaxy 
clusters and groups candidates at z > 1 in the Sp itzer In- 
frare d Camera (IRAC) Shallow Survey (Eiscnhar dt et al] 
I2OO4I ). 

Obtaining spectroscopic measurements for a large num- 
ber of faint cluster galaxies at z ~ 1 is a very difficult and 
time-expensive task. As a consequence, to date there are 
no reliable LFs measured at 2 ~ 1 using only spectroscopi- 
cally confirmed members. Several authors have attempted to 
solve this issue by using photometric redshi fts to define clus- 
ter membership an d to compute the LF (jToft et al ] |2004l . 
iTanaka et al.l [20071 ). but their results have been criticised 
due to the unknown spec tral energy distribution (SEP) o f 
galaxies at high redshifts l|Andreon et al.ll2005l . rShethll2007l ). 

Traditionally, the cluster LF is computed applying a 
background subtraction method, which consists of com- 
puting the difference between galaxy counts in the clus- 
ter and control field directions. Several authors hav e devel- 
oped and applied their own subtraction metho ds (jOemlerl 
1 1974 iPimbblet et"aLll2002l . lAndreon et al]|2005l ). However, 
it is difficult to assess whether this method is free of im- 
portant systematic effects due to possible problems in the 
background subtraction method, which needs to select an 
appropriate area, either around the clusters under study or 
from an independent large area photometric sample. 

It has been claimed that infrared luminosities are bet- 
ter suited to measure the LF than their optical counter- 
parts, since the former are comparatively insensitive to the 
star formation history while reflecting the total stella r mass 
ijCharlotl Il996l . iMadau. Pozzetti. fc D ickinson' '2004"). Fur- 
thermore, infrared luminosities have a negligible extinction 
correction and their K-corre ction depends only we akly on 
Hubble type (see for instance iMannucci et al.ll2001^ . 

The main goal of this work is to study the reliability of 
the statistical background subtraction method to recover the 
underlying observer-frame Ks-band LF. Section 2 explains 
the mock catalogues in which we base our analysis. Section 
3 discusses the background subtraction method and shows 



the resulting LF estimates, and Section 4 shows a compar- 
ison between the underlying and recovered LF, and several 
physical quantities inferred from the LF. Our conclusions 
are summarized in Section 5. 

Throughout this paper we assume Qm = 0.3, Qa ~ 0.7 
and Ho = 70 km Mpc'^ . All magnitudes are in the Vega 
system. 



2 DATA 

This work concentrates in particular, on the possibility to 
study the LF of z ~ 1 clusters, in such a way so as to ob- 
tain the accuracy needed to study different properties of the 
galaxy population i ncluding the relative i mportance of dwarf 
vs. gia nt galaxies (iDe Lucia et al.l |2007| . iGilbank fc Baloghl 
l2008al ). or the expected formation redshift of z ~ 1 clus- 
ter galaxies. To date, the most suited cluster survey for 
such a study is the Red sequence Cluster Survey (RCS-1, 
iGladders etahl |2005| ) However, there is only optical infor- 
mation for these clusters, so a study of the LF over the 
infrared wavelength range would require an observing cam- 
paign to obtain these data. The authors of this work are at 
the moment in the process of obtaining such information, so 
that this study will be possible in the near future. 



2.1 Mock catalogues 

In order to test for possible systematic errors afi'ect- 
ing a statistical measurement of the cluster LF, we ap- 
ply this analysis to mock RCS catalogues. We con- 
struct m ock RCS catalogues using the Millennium sim- 
ulation l|Springel. Frenk fc Whit3 I2OO6I ) populated with 
semi -analytic gala xies from the GALFOR M m odel by 



Baugh et al] (|2005l ) (see also ICole et ahllioOO. Bower et al 



2006; for an alternative approach see lLagos. Cora fc Padilla 



2008^ . The numerical simulation follows the gravitational in- 
teraction of 10" dark-matter particles in a ACDM Universe 
characterized by the cosmological parameters Q,m ~ 0.30, 
Qa = 0.70, a power law spectral index n = 0.94, an ampli- 
tude of matter density fluctuations in spheres of 8 h~^ Mpc, 
as — 0.8, and a Hubble constant Ho = 100 /ikms~^ Mpc~^ 
with h — 0.70. The GALFORM semianalytic model was ap- 
plied to the merger trees of halos in the simulation, which 
have a minimum mass of 1.2 x 10^" h~^ Mq for a minimum 
of 20 dark-matter particles (courtesy of the Durham group). 
The galaxies in the simulation are characterized by their lu- 
minosity in several photometric bands, including the near 
infrared Ks-band and those available in the RCS dataset. 
In order to construct each mock catalogue, we place an ob- 
server in the simulation, and record all angular positions of 
galaxies down to the RCS magnitude limit within an area 
of 100 square degrees. The intrinsic completeness limit in 
our mock catalogues is Ks = 22.0, observer-frame, following 
the resolution constraints of the semi-analytic galaxy pop- 
ulation. Notice that this limit can be lowered to meet the 
restrictions of a given observational sample. 
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Figure 1. Mass function of simulated and observed clusters. The 
filled circles show the mass function of simulated clusters selected 
from the mock catalog. The open triangles show the mass func- 
tion of observed clust ers selected from the RCS-1 cluster catalog 
dGladders et al.ll2007l) . 



In order to include the evolution of the galaxy popula- 
tion with redshift, we use simulation outputs corresponding 
to four different redshifts from z = 0.3 to z = 3 for con- 
secutive redshift ranges. This ensures that the population of 
background galaxies behind our mock RCS clusters are char- 
acterized by an evolving luminosity function to distances 
beyond the maximum expected galaxy redshifts allowed by 
the magnitude limit cut. The fact that we are not interested 
in the differential evolution of the luminosity function with 
redshift justifies our choice of mock light-cones with discre- 
tised evolution. 

The advantage of using mock catalogues relies in that 
the underlying population of clusters is known, along with 
the LF parameters for the field and cluster environments for 
each individual mock. Furthermore, as we are using multi- 
ple mocks extracted from practically independent volumes 
in the simulation, we can estimate the expected effects of 
sample variance in our statistics. Besides these advantages, 
we are also able to make a comparison between the prop- 
erties of galaxies populating the mocks and those in real 
datasets published in previous works, providing a new test 
for the GALFORM semi-analytic model. 

Since we are interested in employing these mock cata- 
logues to study the accuracy of the background subtraction 
method for computing the galaxy LF for RCS-1 clusters, 
we compare the mass function between simulated and ob- 
served clusters. Figure [1] shows the mass function of simu- 
lated clusters selected from t he mock catalogues and RCS-1 
observed clusters taken from iGladders etahl (I2007D . where 
it can be seen that both mass functions are roughly com- 
patible for M200 > 1-5 X 10^* Mq/i"! to a lo- level. The 
offset in the observational data can be due in part to errors 
in the mass measurements wh ich usually broaden statistical 
distributions (see for instance iPadilla fc Lambai (|l999l )) 



3 LUMINOSITY FUNCTION ESTIMATES 

The optimal method to build a reliable cluster LF con- 
sists in using spectroscopic information to determine cluster 
membership. However, obtaining spectroscopic redshifts for 
a large number of faint galaxies at high redshifts is pro- 
hibitively expensive in telescope time. 

The cluster samples in our mock catalogues are charac- 
terized by a z = 1 LF in the Ks-band which we fit with a 
Schechter function of characteristic luminosity Ks* = 18.27 
and faint-end slope a = —1.36. Our approach to deter- 
mine the reliability of different LF estimators takes ad- 
vantage of our knowledge of the underlying LF parame- 
ters. We notice that the characteristic luminosity in the 
model is in good agreement with the observa tional result 
Ks* = 18.5 ± 1.2 found by IStrazzullo et all ^QM), who 
applied a background subtraction method and fitted for 
both, the characteristic luminosity Ks* and a, the faint-end 
slope. Other measu rements show d i fferen t levels of agree- 
ment; for instance, iKodama et al] (|2003D found a fainter 
value Ks* — 17.57 ± 0.24, using a fixed value of a = —0.9 
and relatively shallow nIR data. With respect to the faint- 
end slope, this parameter is poorly constrained by observa- 
tional data, and consistent with a = — 1. Therefore, down to 
the accuracy of the observational results, the semi-analytic 
model used to construct the mock catalogues is able to re- 
produce the observed cluster galaxy population at 2 ~ 1. 

In the following subsection we discuss two possible LF 
estimators that we will apply to the mock cluster catalogues. 



3.1 Background subtraction and photometric 
redshift LF estimators 

Two widely used methods to calculate the LF are the statis- 
tical background subtraction and the photometric redshift 
method. The former consists on using one or several con- 
trol fields to determine the number of contaminating galax- 
ies (background and foreground galaxies) per unit area as 
a function of magnitude, to then compute the galaxy LF 
as the difference between the galaxy number count in the 
cluster direction and the control field. This is the method of 
choice when there is photometry in only one or two bands for 
galaxies in the cluster direction, and very few cluster mem- 
bers have been spectroscopically confirmed. Employing nu- 
merical simulations, it has been shown that the background 
subtraction method can recover accurately t he underlying 
LF o f clusters selected in three dimensions (jValotto et al.l 
120011 ). 

The background subtraction method requires to define 
a magnitude bin width, a cluster-centered aperture, and a 
control field region. When computing the number counts in 
the cluster direction for the RCS mock catalogue, we mimic 
the observational method of employing galaxies enclosed by 
a circular aperture of radius 0.27 h~^ Mpc centered on 
the brightest cluster galaxy and a magnitude bin of width 
0.5 mag. This aperture size is justified because the nIR imag- 
ing of RCS clusters covers only 0.8 h~^ Mpc of the cluster 
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central region. We defined the control field region based in 
the observational nIR surveys to date, which should have a 
minimum 5a limiting magnitude of Ks = 22.0 and cover a 
big area in the sky. The candid ates were the Faint Infrared 
Extragalactic Survey (FIRES: ITabbe et al.l 120031 '). which 
covers 8.3 arcmin^ and reaches a 5a depth for point sources 
of Kssct = 24.5; the Great Observatori es Origins Deep Sur - 
vey of CDF-S region (GOODS /IS AAC: lRetzlaff et al.ll2008h , 
which covers 170 arcmin'^ and reaches Ksso- = 22.8; and the 
nIR VIMOS-VLT Deep Survey (VVDS-NIR; llovino et all 
I2OO5I ). which covers 430 arcmin and reaches Ksso- = 22.1. 
Figure [5] shows the ratio between the rms variation of num- 
ber count along the control field regions previously men- 
tioned and the rms variation of actual cluster members in 
the mock. We can see that a control field of size the FIRES 
survey contributes about 30% to the rms variation along 
the cluster direction, so the accuracy of the LF built from 
a background subtraction method would be highly limited 
by the small control field region. Also, in order to reach a 
0.4% contribution of background rms variation to the Ks- 
band LF error bars we would need a nIR survey covering 
1 deg^, which does not even exist at this magnitude depth 
to date. This way, GOODS/ISAAC and VVDS-NIR are the 
most suitable surveys to build the Ks-band LF with about 
100% completeness to magnitude 21.0 using the background 
subtraction method. 

We discard the VVDS-NIR survey, since the calibrated 
images of GOODS are publicly available and are taken with 
the same instrumen t as the nIR imaging of the RCS clusters 
l|Mufioz et al.ll2008l ) 

Using the mock catalogues, we compute the rms vari- 
ation of number counts from several 170 arcmin'^ regions 
in the mock catalogues, and compare this to Poisson errors 
and jacknife algorithm estimates (see Figure [S] top panel). 
It can be seen that a Poisson distribution underestimates 
the actual rms variation about 30% and 60% for bright and 
faint galaxies, respectively. This high discrepancy can be ex- 
plained by the galaxy-galaxy correlation and the presence of 
large sc ale structure along the control field region (see for in - 
stance, iPadilla et al]|2004l . |Paz. Stasvszvn fc Padiliall2008h . 
We also see that the jacknife algorithm reproduces the un- 
derlying spread in number counts seen in the simulation to 
a higher degree. This way, the error in the galaxy count 
along the control field region, ajv^:, will be computed from 
the jacknife algorithm. 

The observational program being c arried out to obtain 
near-infrared imaging of RCS clusters (|MufLOz et al.|[2008l ) 
will only cover fields centered in clusters, therefore, it is im- 
portant to test whether the use of number counts from a 
different source catalogue will still allow measuring the LF 
with minimum systematic deviations from the underlying 
values. 

The number of clusters members along the cluster di- 
rection in the j*'' magnitude bin Nj , and its error, crjv , are 
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Figure 2. Ratio between the rms variation of number counts 
along the control field regions of several sizes and the rms vari- 
ation of actual cluster members. Open circles correspond to the 
deepest nIR survey FIRES, stars correspond to the public ESO 
survey of the CDF-S region GOODS/ISAAC, and triangles to the 
deep and wide VVDS. The squares correspond to an hypotheti- 
cal nIR survey covering 1 deg^, reaching a 5a limiting magnitude 
deeper than Ks = 22.0 



given by 

A^, = Ndj - Nfj 



aNi 




(1) 
(2) 



where Ndj is the number of galaxies (cluster-|-field) in the 
cluster direction in the j*'' bin, Nfj is the number of galaxies 
in the control field region in the j"" bin, Qd is the area of 
the cluster-centered circular aperture, and f2/ is the area of 
the control field region. The error in the number counts in 
the cluster direction, o-jv^^ , is estimated by the Poissonian 
error, a choice we justify next. 

We calculate the counts in the cluster direction for 100 
independent samples containing 10 RCS clusters each, and 
calculate the dispersion in these counts as a function of 
apparent magnitude. We then select one of these samples 
at random, which we use to compute the classical Pois- 
son error \/iV for the number counts in the cluster direc- 



tion, and the mod if ied P oisson error 1 -|- a/ {Ndj + 0.75) pro- 
posed by GehreQ (Il981) and employed by s everal authors 
lAndreon et all |2005| . IStrazzullo et al.l 120061 ). The bottom 
panel of Figure |3] shows the comparison between the disper- 
sion in the counts in the cluster direction and its estimated 
error as a function of the observed magnitude. From this 
figure we conclude that the classical Poisson error resem- 
bles best the dispersion in the mocks, while the modified 
Poisson errors overestimates errors by a significant amount. 
Therefore, this justifies our choice of Poisson errors when 
study ing the cluster counts in the mock CLFs. IZheng et al.l 
( 12005! ) find that the number of galaxies in dark matter haloes 
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Figure 3. rms variation in the galaxy counts along the control 
field region of size the GOODS/ISAAC survey (top panel) and 
along the cluster direction (bottom panel) as function of magni- 
tude. The solid line shows the actual dispersions built from 100 
independent mock samples. The dashed and dotted lines in the 
top panel show the errors in the background number count esti- 
mated from the classical Poissonian error \/N and from applying 
a jacknife technique, respectively. The dashed and dotted lines 
in the bottom panel show the errors in cluster direction num- 
ber count estimated from the classical and a modified Poissonian 
error 1 -|- V^V+TX75 



in semi-analytic and SPH simulations follows a Poisson dis- 
tribution, which would indicate that galaxies outside the 
cluster (either in front or behind due to projection effects) 
only contribute a relatively low fraction of the uncertainty 
in the cluster counts. 

The photometric redshift method consists of using 
broad-band photometry in several band-passes as a very low 
resolution spectrum, and then fitting a spectral energy dis- 
tribution (SED) template to estimate the redshift of the 
galaxies in the cluster field. The final step consists on defin- 
ing as cluster members those galaxies lying within a red- 
shift shell centered on the cluster redshift. For simplicity 
, we do not compute the photometric redshifts of 
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Figure 4. Completeness {solid line) and failure ratio [dashed 
line) of photometric redshift method as function of Ks-band mag- 
nitude. The open circles correspond to cluster members candi- 
dates lying 1 crz/{l+z) from the cluster redshift, while the trian- 
gles and squares correspond to galaxies lying 2 and 3 (^2/(1+2) > 
respectively. 



galaxies from fitting a SED to the observed magni- 
tudes but rather from perturbing their actual red- 
shifts by a gaussian distribution with ^^/(i-i-z) — 0.05 . 
It is important to note that the true probability dis- 
tributions of photometric redshifts are often highly 
non-Gaussian, and that secondary maxima or ex- 
tended wings might change the accuracy of photo- 
metric redshift method used in this paper. Figure |3] 
shows the completeness and failure ratio of the photomet- 
ric redshift method computed from 500 independent Monte- 
Carlo realizations of photometric redshift of randomly se- 
lected clusters. Each realization consists in choosing 
a random cluster located at 2 ~ 1, then selecting 
those galaxies within a circular aperture of radius 
0.27 h^^ Mpc, and finally, simulate their photometric 
redshifts using the actual redshifts from the mock 
catalogues via a Monte-Carlo procedure. We can see 
that the photometric redshift cut l-cr2/{i_|_2) recovers a low 
percentage of actual cluster galaxies and has a failure ra- 
tio of about 10 percent, while the redshift cut of 2-a2/(i+z) 
recovers almost all the actual cluster galaxies and has a fail- 
ure ratio of about 15 percent. The cut of 3-crj./{i_|.2) has the 
highest completeness, but the failure ratio of bright galaxies 
increases rapidly. According to the previous considerations, 
we define as cluster members those galaxies lying within 
2-o"z/(i+z) from the cluster redshift. Large scale structure 
plays a role in this calculations and therefore a simple ana- 
lytic gaussian is only approximate. 
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The RCS clusters under study by iMunoz et all (120081 ) 
have deep J and Ks-band imaging, and shallow Rc and z 
imaging, so it is not possible to build the LF up to mag- 
nitude Ks = 21.0 from applying the photometric redshift 
method. In order to take advantage of the available data, 
we propose a new method that consists in using photomet- 
ric redshift of galaxies to compute the bright part of the 
LF and using the background subtraction method to com- 
pute the faint part. The cluster members with Ks ^ 18.5 
are identified by their photometric redshift, while the num- 
ber of members with Ks > 18.5, by using the background 
subtraction method. The combination of both methods will 
be called B+Z method hereafter. 



Table 1. Best-fitting parameters for Ks* and a as free pa- 
rameters. 



Dataset 


Par am 


Actual 


Background 
method 


B+Z 
method 


CSl 


K* 
a 


18.27 
-1.36 


18.95«:^» 


18 94+0-''6 

-i.io«:2 


CS2 


K' 
a 


18.27 
-1.36 


17.98+0-*? 

_1 45+0.30 


18.56l°:« 
-1 29+0-28 



CSS K* 18.27 17.71+0;;5' I8.I8IO;: 

a -1.36 -l-50+g:t2 "l-^^^o. 



3.2 The composite luminosity function 

In most galaxy clusters there are too few galaxies to deter- 
mine accurately the shape of the LF. A solution to this prob- 
lem is to measure the LF using a combined sample of clus- 
ters, allowing a higher precision measurement. This method 
is commonly known as the Composite Luminosity Function 
(CLF) method. 

The CLF is built according to the following formulation. 



1 



E 



(3) 



where Ncj is the number of galaxies in the j-th bin of the 
CLF, Nij is the number of galaxies in the j-th bin of the 
i-th cluster, and rij is the number of clusters with limiting 
magnitude deeper than the j-th bin. The formal errors of 
the CLF were computed according to, 



1 



clustet 

E 



(4) 



where ajv^j is the formal error of the number of galaxies in 
the j-th bin of the i-th cluster. 

For computing the CLF, we selected several galaxy clus- 
ters with masses about 1.5 x 10^** M0, located between red- 
shifts 0.90 < z < 1.10. The observed CLF was built for 3 
cluster samples taken from the mocks: the first sample con- 
tains 5 clusters (CSl), the second, 10 clusters (CS2), and the 
third, 50 clusters (CSS). We reset the completeness limit in 
our mock catalogu es to Ks = 21.0 in order to match the 
RCS nIR imaging (|Munoz et al.ll20o"^ '). and we apply both, 
the background subtraction and B-l-Z methods for comput- 
ing the CLF. 

Figure [5] shows the resulting CLF computed for a CSS 
dataset, using both the background subtraction and B-l-Z 
methods (upper and lower panels, respectively). The open 
symbols correspond to the underlying CLF calculated using 
all the cluster members, and the filled circles correspond to 
the CLF obtained using only data available in observational 
samples. This figure will be discussed into more detail in the 
following section. 



Table 2 . Best-fitting parameters derived with a fixed value of 
a = -1.36 



Dataset 


Par am 


Actual 


Background 
method 


B+Z 
method 


CSl 


K* 


18.27 




18.50lo2g 


CS2 


K' 


18.27 


18.30«:- 


18.4210:1^ 


CSS 


K* 


18.27 


18.20«:0^ 





4 ANALYSIS 

4.1 Schechter parameters 

In order to study the reliability of the background subtrac- 
tion and B+Z methods to recover the underlying LF, we fit 
a Schechter function to the observed CLF and then com- 
pare the best fitting observed parameters to those of the 
underlying CLF. 

The best fitting parameters are found using a minimum 



X method given by. 



E 



N{mj 



where N(mj 



(f>(m)dm, 



(5) 



(6) 



; a) is the predicted value of the CLF at 
the magnitude bin rrij, mi,i„ is the magnitude bin size, and 
(jjim) is the Schechter function. 

We search the global minimum in the full (j)* — ra* — a 
space sampled by a parameter grid. The best-fitting param- 
eters after marginalizing over </>* , for Ks* and a taken as free 
parameters, are summarized in the Table[T]for the three sam- 
ples of mock galaxy clusters defined in the previous Section. 
The systematic errors are defined as the difference between 
the actual and estimated value, while the stochastic errors 
are defined by the 68 percent confidence levels. 

Figure |S] shows the 1- and 2-a likelihood contours of 
the Schechter function parameters for our sample CSS us- 
ing both, the background subtraction and B+Z methods. 
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Figure 5. Ks-band composite luminosity function (CLF) built from the background subtraction and B+Z methods. Shown are the CLE 
built from the background subtraction method for 5 clusters (top-left), for 10 clusters (top-right), and for 50 clusters (bottom-left) . Also, 
the CLF built from the B+Z method for 50 clusters (bottom- right) . Dashed line shows the underlying CLF, while open squares show the 
recovered CLF from applying the respective method. 



The plus symbols mark the actual CLF parameters and 
the crosses mark the best fitting observed CLF parameters. 
From this figure, we can conclude that the B+Z method 
gives more accurate and better constrained results than the 
background subtraction method, since the stochastic error 
of Ks* is reduced by almost a factor or 2. 

It is useful to study the reliability of the best LF mea- 
surement method (B+Z estimator) as a function of the num- 
ber of galaxy clusters used for computing the CLF, since this 
can have a direct impact on current and future observational 
campaigns. In principle, a larger sample of clusters should 
put better constraints on the values of Ks* and a as well as 
lower systematic errors. Figure [7] shows the stochastic and 
systematic errors for both Ks* and a. From this figure we 
can conclude that the method reliability and accuracy de- 
pend strongly on the number of galaxy clusters, as increasing 
the sample size by a factor of 10 reduces the stochastic and 
systematic errors in Ks* by factors of 3 and 30, respectively. 

When the parameter a is fixed to the true underlying 
value, errors also improve, for instance, by a factor of 3 in 



the characteristic luminosity. The results for the three sam- 
ples of mock galaxy clusters are summarized in Table [S] In 
the case where the value of a is unknown and a fiducial value 
a = — 1 is assumed (as may be the case in an observational 
estimate of the CLF), the systematic errors increase signifi- 
cantly due to the likelihood degeneracy in the Ks*-a plane, 
which in our mocks forces the resulting characteristic lumi- 
nosity to a fainter value. It is also noticeable that in this case, 
the stochastic error is lower than when fixing a to the true 
underlying values; this indicates that even though the un- 
certainties in Ks* diminish when fixing a, if this parameter 
is offset by only a 23% of its value, there could be significant 
systematic errors in Ks* producing biased analyses. 

Regarding the faint-end slope, a, the stochastic and sys- 
tematic errors also diminish by about a factor of 4 when in- 
creasing the sample from 5 to 50 clusters. It is important to 
notice that regardless of the number of clusters in the sam- 
ple, the stochastic error is always found to be larger than 
systematic deviations. 
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Figure 6. Contour plots showing tlie constraints on the Schechter 
function parameters derived from the minimization analysis 
applied to the CLF obtained using 5 clusters {dashed line) and 
50 clusters {solid line). The upper panel shows the results from 
applying the statistical background subtraction method; the lower 
panel shows the results from applying the B+Z method. The plus 
sign marks the parameters of the true underlying CLF, and the 
cross marks the best fitting parameters of the measured CLF. 
The labeled contours represent the 1- and 2-a confidence levels. 



4.2 Mean formation epoch 

The study of the evolution of the CLF is a powerful tool to 
improve our knowledge on the processes of galaxy formation 
and evolution, in particular in high density regions of the 
Universe where the probability of galaxy interactions with 
their environment is at its highest. In addition, these studies 
are easier to interpret when using infrared luminosities, since 
these are better suited than optical photometry to measure 
the LF, since the former reflect the total stellar mass of 
galaxies and do not depend strongly on the details of their 
stellar populations (jGavazzi. Pierini fc Bosellilll996l ). 

One of the most interesting quantities that can be in- 
ferred from the evolution of the LF is the mean redshift of 
the bulk of star formation in cluster galaxies, zf, this red- 
shift is defined as a measurement of the time when these 
galaxies have already acquired most of their stellar mass. 



Figure 7. Upper panel: Stochastic and systematic errors in Ks* 
obtained through the B+Z LF measurement method, for samples 
composed by difi^erent number of galaxy clusters, for a fixed or free 
value of the a parameter (line types are shown in the figure key). 
Upper and lower dotted lines show the stochastic errors in Ks* 
for fixed values of a = —1.36 and a = —1.00 respectively. Upper 
and lower dot-dashed lines are the systematic error in Ks* for 
a = —1.00 and —1.36, respectively. The horizontal long-dashed 
line is the maximum error that still allows to distinguish between 
passive evolution scenarios with formation redshifts of z j = 2 and 
Zf = 5. Lower panel: Stochastic and systematic errors in a. 



This is most easily obtained from observational estimates 
of the evolution of Ks* by direct comparison to predictions 
using passive evolution models with different Zf. 

We employ th e pa ssive evolution models of 
iKodama fc Arimotol ()l997l ) to compute the minimum 
number of galaxy clusters needed to distinguish between 
a scenario with Zf — 2 and Zf = 5 for the galaxies in 
our mock catalogues. According to these models, for a L* 
galaxy observed ai z — 1, the difference in Ks* between 
both scenarios is about 0.3 mag. Therefore, in order to 
compute Zf with an accuracy of 2 Gyr at the 68 percent 
confidence level, the maximum systematic and stochastic 
errors allowed on Ks* would be AKs* — 0.15 (shown as a 
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Figure 8. Dwarf-to-giant ratio (DGR) for red-sequence cluster 
galaxies computed from background subtraction {solid line) and 
photometric redshift {dashed line) methods. The triangles and 
squares correspond t o the DGR est i mated from applying the 
magnit ude limits by iDe Lucia et alj ll2007h and ICilbank et al.l 
ll2008lj) , respectively. The underlying value is shown by the dotted 
line. The error-bars sho w the stochastic errors. A small horizontal 
offset was applied to the lGilbank et al.l ll2008d) results to improve 
clarity. 

horizontal long-dashed line in figure [T]! . According to the 
results showed in figure[7]for the B+Z method, leaving both 
Ks* and a as free parameters, this would require a sample 
with a total number of 520 clusters at 2 ~ 1. This number 
can be lowered by fixing a to its actual value —1.36, as 
shown by the upper dotted line in figure [71 since this would 
then require only 10 clusters in total. However, fixing a to 
the flatter value a — —1.00, as shown by the lower dotted 
line in figure [T] would not allow to compute an accurate 
value of 2/ and a comparison between formation scenarios 
would be biased. Another option could be to combine 
the measurement of Ks* from an available observational 
sample, with alternative methods to measure the formation 
epoch of galaxies such as Bruzual & Chariot (2003) fits to 
the SEDs of confirmed members oi z = 1 galaxy clusters, 
with possible restrictions from photometric data depth and 
number of galaxies with such information. 



4.3 Dwarf-to-Giant ratio 

The dwarf-to-giant ratio (DGR) of red-sequence galaxies 
is a commonly used parameter that allows an interpreta- 
tion of the cluster galaxy population at a given redshift, 
providing the relative importance of the dwarf and giant 
galaxy population in clusters. This parameter is simply de- 
fined as the ratio between the number of faint and bright 
red-sequence cluster galaxies within predefined magnitude 
ranges. As such, it does not require an analytic fit to the 
measured LF and therefore does not need forcing any par- 
ticular function to the real distribution of luminosities, while 
also avoiding dealing with parameter degeneracies. However, 



it does require to define arbitrary luminosity limits for the 
dwarf and giant populations that usually depend on the 
available data. A good revie w about cosmic ev olution of the 
DGR was published by Gilbank et all l|2008ah . 

Since we are interested in studying the accuracy and 
bias of the DGR value computed from applying the back- 
ground subtraction method, we compute the actual and esti- 
mated DGR values from the mo ck catalogue. W e employ the 
magnitude limits published by iDe Lucia et al] (|2007i ). who 
studied the red-sequence galaxy LF and defined as luminous 
galaxies those satisfying My ^ —20.0, and as faint galaxies, 
-20.0 < s; -18.2 (hereafter DEL07). F or comparison 
we also use the magnitude limits proposed bv lGilbank et al] 
(,2008bl ). who defined luminous galaxies as those within the 
range —22.7 < Mv ^ —20.7, and faint galaxies those 
within -20.7 < My ^ -19.7 (hereafter GIL08). We con- 
vert these rest-frame V-band magnitudes t o our observer- 
frame Ks-band, using the same method as iDe Lucia et al.l 
l|2007l ). We use the GALAXEV stellar population synthesis 
code (|Bruzual fc Chariot! lioo^ ') to generate model galaxy 
SEDs arising from a single burst stellar population formed 
at Zf = 3. 

In order to compute the actual value of the DGR, we 
build the CLF for cluster g alaxies in the red-se quence from 
100 mock clusters. In the iBaugh et al.l (|2005l ) model, the 
z = 1 rest-frame color-magnitude relation (CMR) is defined 
by the mean color (r — z)rest = 2.89 — 0.025 x Zrcst and 
intrinsic dispersion A(r — z) = 0.13. We obtain values of 
DGRdel = 2.94 and DGRgil = 1.83 for the magnitude 
limits defined by DEL07 and GIL08, respectively. 

The observational estimates found by DEL07 and 
GIL08 correspond to DGR = 1.05 ± 0.16 ai z ~ 0.75 and 
DGR — 0.84 ± 0.05 at z ~ 0. 9, respectively; as can be 
seen, the semi-analytic model bv lBaugh et al.l l|2005h tends 
to produce a larger number of dwarf galaxies at a redshift 
of 2; ~ 1, and given the observed evolution of DGR, the dis- 
crepancy between the theory and observations will be even 
higher at this redshift. Using the paramete rization of the 
evolution of the observed DGR presented bv lGilbank et all 
(:2008ai), the expected value of DGR at z = 1 for the magni- 
tude limits defined by DEL07 is about 0.88, which is a factor 
3.3 smaller than the value computed from our mock cata- 
logues. Notice that the value of the faint-end slope in the 
model, a = —1.36, is probably responsible for this discrep- 
ancy, since observed values tend to favour a > — 1, although 
with large uncertainties. 

We use the underlying DGR values to study the sys- 
tematic and stochastic errors in the DGR as a function 
of the number of clusters in the samples studied. Figure [S] 
shows these results, where triangles and squares were com- 
puted using the DEL07 and GIL08 magnitude limits, re- 
spectively. The solid lines correspond to the DGR computed 
from the background subtraction method, while the dashed 
lines were computed from the photometric redshift method 
as explained in section 3.1 It can be seen that systematic 
errors in both DEL07 and GIL08 are small compared to the 
dispersion in the measurements, and also, that stochastic 
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Table 3. Cluster occupation numbers derived from the mea- 
sured CLFs. 



Dataset 


Measured value 


Actual value 


CSl 


0.80 ±0.42 


0.61 


CS2 


0.73 ±0.30 


0.61 


CSS 


0.76 ± 0.07 


0.61 



errors show a clear improvement as the number of clusters 
increases. For the background subtraction method, the lat- 
ter goes from about a 35% percent error for samples of 5 
clusters to about a 5% percent error for samples with 50 
clusters. 



4.4 Occupation numbers 

A commonly used quantity in the study of galaxy popu- 
lations in clusters is the occupation number, or the total 
number of galaxies above a lower limit in luminosity in each 
galaxy cluster of a given mass. This can obtained from the 
analysis of galaxy luminosity and correlation functi ons in a 
form a lism called the Halo M odel (see for instance ICooravl 
I2OO5I . ICoorav fc ShethI \200±. A study o f the correlation 
function of SDSS galaxies by Zehavi et al.l (|2004) shows that 
« ~ clusters of galaxies of mass 1.5 X 10^*h-^Mpc would 
contain an average of ~ 3 galaxies brighter t han Mr = — 21, 
equivalent to M* — 0.5 l|Blanton et aLllioOSl ). The same cal- 
culation performed on the z = galaxy population in the 
iBaugh et al] (|2005l ) model indicates that clusters of the same 
mass in the model are populated by an average of 2.2 ±0.1 
galaxies, in reasonable agreement with the observational es- 
timates. 

Estimates of the LF of clusters of a given mass can 
be used to estimate this number at any given redshift. In 
our study we use the estimated counts down to Ks* — 0.5 
to find the number of galaxies residing in z ~ 1, simu- 
lated RCS clusters. The resulting numbers depend on the 
number of clusters available in the sample as is shown in 
Table [3] As can be seen, the measured occupation num- 
bers are in agreement with the underlying values for CSl 
and CS2 samples, although the measurement errors dimin- 
ish significantly as the number of clusters in the sample in- 
creases. The underlying value of 0.61 galaxies per cluster 
does not change significantly when calculating this quantity 
using a different photometric band, since the lower lumi- 
nosity limit depends on the characteristic luminosity of the 
CLF; therefore, the occupation number of simulated clusters 
of M ~ 1.5 X 10^ 
z — 1 to z = 0. 



'h ^M0, increases by a factor of ~ 3.6 from 



5 CONCLUSIONS 

We studied the reliability of different LF measurement meth- 
ods appropriate for the study of high redshift, z ~ 1, clus- 
ters of galaxies. We do this by using mock cluster cata- 



logues constructed using GALFORM (|Baugh et all |2005| ) 
semi-analytic galaxies that populate a numerical simulation 
of a ACDM Cosmology. Our mock catalogues mimic the RCS 
clusters and the GOODS/ISAAC survey, which are the tar- 
gets of an ongoing observational project that will comple- 
ment RCS cluster photometry with Ks-band data, and will 
use the GOODS/ISAAC survey t o obtain the backgr ound 
contamination in cluster regions (jMufioz et al. liooi). We 
find that the best way to characterize the error in the back- 
ground counts is by means of the jacknife method for at least 
4 jacknife subsamples, which takes into account the varying 
structures in the line of sight of pencil-beam surveys. 

Our studies indicate that the joint use of photomet- 
ric redshifts of bright galaxies plus background-corrected 
counts of faint galaxies (B±Z method), provides the best 
recovery of the underlying CLF. In the process we also 
found that the classical Poisson method provides accurate 
errors for the galaxy counts in the cluster direction. The er- 
ror in the best-fitting Schechter parameter Ks* from B±Z 
method, decreases by almost a factor of 2 when compared 
to the measurements obtained from using a method relying 
on background-corrected counts alone. 

Even though the B±Z method provides the best re- 
sults, the measured CLF shows large stochastic errors for 
datasets with small numbers of galaxy clusters. For a dataset 
of 5 galaxy clusters we find 4% and 44% stochastic errors 
in Ks* and a, respectively (AKs* ~ 0.7 magnitudes, and 
Aa ~ 0.26). However, increasing the sample size to 50 clus- 
ters improves this result in a dramatic way, reducing these 
errors to a 2% and ~ 9% (AKs* ~ 0.1 magnitudes, and 
Aa ~ 0.05). The use of a fixed value of the faint-end slope 
a for the Schechter fit improves the stochastic errors by a 
factor of 3; systematic errors do not change significantly, 
unless the assumed value of a is offset from the underlying 
value. 

The accuracy of the B±Z method in recovering the un- 
derlying dwarf-to-giant ratio (DGR) depends strongly on the 
luminosity limits used to define the bright and faint galaxy 
population, but in general the systematic error does not 
improve significantly when using larger samples. However, 
there is a clear tendency to obtain smaller stochastic errors 
in this quantity. On the other hand, in order to distinguish 
between formation redshifts of Zf — 2 and 5 (an uncertainty 
of ~ 2Gyr) for the bright cluster members, the sample of 
clusters needs to contain a minimum of ~ 10 clusters; this 
result corresponds to the B±Z method for measuring the 
CLF, plus leaving only the characteristic luminosity Ks* as 
a free parameter and fixing a to an appropriate value. The 
occupation number of clusters in the sample, that can be 
associated to the median cluster mass, can be accurately 
obtained when using 10 or more clusters in total. The ac- 
curacy of this measurement increases with the sample size 
reaching a stochastic error ~ 3% for a total of 50 clusters. 
Such a measurement can also be used to study the evolu- 
tion of the halo occupation number as a function of redshift, 
which can help constrain galaxy evolution models. 

In conclusion, the clear advantage of near infrared wave- 
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lengths for the study of galaxy luminosities in clusters at 
relatively 2: ~ 1 can be exploited to a maximum if, i) multi- 
band photometry allowing a photometric redshift estimate 
for the brightest cluster members is available, so that a joint 
method using background subtraction and photometric red- 
shifts can be applied to calculate the CLF, ii) a reasonable 
number of > 50 clusters is included in the sample, and iii) 
the available measurements of CLFs at these redshifts al- 
lows to assume a fixed value for the faint-end slope of the 
luminosity function a. 

Regarding the modeling of the evolution of galaxies in 
the semi-analytic model, the current precision in the CLF 
measurements does not allow us to reach firm conclusions 
on whether the model succeeds to reproduce the observed 
galaxy population. There is, though, an indication of an 
excess of dwarf galaxies by approximately a factor of 3 in 
the model with respect to observations. However, we have 
demonstrated that only a factor of ~ 10 increase in sample 
size with respect to those used in currently available Ks- 
band CLF measurements, would make a more stringent and 
decisive test possible, allowing the use of cluster galaxies at 
2 ~ 1 to further our understanding of how galaxies form and 
evolve in the Universe. 
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